Ab initio correlation approach to a ferric wheel-like molecular cluster 
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We present an ab initio study of electronic correlation effects in a molecular cluster derived from 
the hexanuclear ferric wheel [LiFe6(OCH3)i2-(dbm)6]PF6. The electronic and magnetic properties 
of this cluster have been studied with all-electron Hartree-Fock, full-potential density functional cal- 
culations and multi-reference second-order perturbation theory. For different levels of correlation, a 
detailed study of the impact of the electronic correlation on the exchange parameter was feasible. 
As the main result, we found that the influence of the bridge oxygen atoms on the exchange pa- 
rameter is less intense than the influence of the apical ligand groups, which is due to the geometry 
of the cluster. With respect to the cluster model approach, the experimental value of the exchange 
parameter was affirmed. 
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I. INTRODUCTION 

In contemporary condensed-matter physics, the role 
of molecular magnetism is steadily growing and receives 
attention from both experimental and theoretical physi- 
cists and chemists 0, S IE 13- The most studied and 
therefore best known molecule in this field is the Mn- 
12-acetate J5J], but the class of the ferric wheels (wheel- 
shaped iron rings) is becoming an important subject of 
various studiesjy, S 13 ■ 

In the past few years, molecular magnets and some 
ferric wheels have been treated with ab initio quan- 
tum chemical methods like the Hartree-Fock approx- 
imation and density functional approaches using sev- 
eral functionals like the local density approximation 
or hybrid functionals. The usually observed results 
are that Hartree-Fock theory often strongly underesti- 
mates physical properties like the exchange parameter 
[9l Hoi flu Q2I . while the DFT methods overestimate those 

values HIH El E3 

This paper will deal with a molecular cluster de- 
rived from the hexanuclear ferric wheel [LiFeg(OCH3)i2- 
(dbm) 6 ]PF 6 |p,llal- The previous analysis based on the 
full molecule Ij showed that one-determinantal Hartree- 
Fock theory failed to reproduce the observed exchange 
parameter J— -21 K 16]. Density functional calculations, 
on the other hand, indicated that there is an enormous 
dependence of the computed exchange parameter on the 
functional chosen. This problem was already observed 
earlier for other systems, e.g. [l5|, H3, I3- A possible 
solution to this discrepancy is to consider the electronic 
correlation, which is strongly influencing the exchange 
parameter, by wave function-based methods. This should 
lead to a more controlled description of the magnetic be- 
havior of the compleXjjas it was demonstrated for various 
systems (e.g. HI \%L M, H El El E3 ) • 



FIG. 1: Geometry of the molecular cluster. The model 
complex was derived from the full molecule [LiFe6(OCH3)i2- 
(dbm)6]PF6 and was slightly modified to achieve C a (x)- 
symmetry. 
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For antiferromagnetic systems like the ferric wheels, 
the magnetic exchange can be qualitatively described 
considering the electronic charge transfer between the 
magnetic centers over a bridging atom or between the 
magnetic centers and various ligand groups. Within a 
cluster approach, the electronic charge transfer can be 
analyzed within a multi-reference second-order pertur- 
bation theory scheme (MRPT2) [28|. MRPT2 is very 
similar to the complete active space second-order per- 
turbation theory (CASPT2) that has been found to give 
very accurate results for the exchange parameter (see, 
e.g. [H El E3). The MRPT2 method is based on a 
reference ground state wave function which can be cho- 
sen as a multiconfiguration self-consistent field (MCSCF) 
wave function. This reference wave function is the initial 
point for treating the electronic correlation perturbation- 
ally. In some cases, the occurrence of intruder states can 
cause convergence problems |23|, l29l. and the level-shift 
technique proposed by Roos et al [3(| has to be applied. 

In this paper, we analyze the electronic properties 
of the simplified complex displayed in Figure ^ The 
model is derived from the molecule [LiFe 6 (OCH 3 )i 2 - 
(dbm)g]PFg whose structure has been determined by Ab- 



bati et al [T(j . We first apply the Hartree-Fock and a hy- 
brid functional (B3LYP) approach, in order to verify the 
cluster model which is used to represent a fragment of the 
full molecule. Then the effect of electronic correlations 
is studied more detailed, with second order perturbation 
theory at the MRPT2 level. The influence of the level 
shift on the MRPT2 results is determined in detail. 



II. METHOD 

The geometry of the molecular cluster is based on the 
measurements of the primal ferric wheel [LiFeg(OCH3)i2- 
(dbm)e]PF6 by Abbati et al [16j. From these data, a 
complex consisting of two iron atoms and some ligands 
was modeled (see Figure QJ, as this is the maximum 
what can be treated by MRPT2. The iron atoms are 
six-fold coordinated and thus have the proper coordina- 
tion of the Fe ions like in the full molecule. Two point 
charges of +1 were added at the position of the neighbor- 
ing iron atoms in order to restore the charge neutrality 
of the cluster. The CgHs rings were replaced with hydro- 
gen atoms (bonding length 0.93A) as well as the methyl 
groups at the bridge oxygen atoms (0.95A). To achieve 
C s (a;)-symmetry which is necessary to keep the MRPT2 
calculations tractable, the positions of the atoms were 
slightly modified. The full geometry is given in Table |U 

For a proper description of the physical properties of 
the molecular cluster, calculations with the codes CRYS- 
TAL2003 HH IP and MOLPRO2002 [33| were carried 
out. Within the scope of the CRYSTAL calculations, we 
employed the unrestricted Hartree-Fock (UHF) method 
and the hybrid functional B3LYP (a functional with ad- 
mixtures, amongst others, of functionals by Becke, Lee, 
Yang and Parr) . Note that the CRYSTAL code can treat 
systems of any periodicity, so that the molecular cluster 
was treated as a single molecule, i.e. not as a periodic 
system. These calculations are of broken symmetry type 
|34l I35I 1361 137| , as the space symmetry is lowered. The 
state is not an eigenfunction of S 2 , but only of S z . 

For the needs of a molecular system, we first chose the 
same basis set as used in ref. [SJ] for the full molecule, 
which was found to be reliable. Henceforth this basis 
set is referred as basis set A. In contrast to the code 
MOLPRO, the possibility of adding point charges is not 
implemented in the CRYSTAL code. A point charge of 
+1 can however be achieved by a H atom with a basis 
function with a very high exponent (100000 a.u.), which 
does not allow charge transfer to the H atom and thus 
acts like a point charge. The idea of the CRYSTAL calcu- 
lations was to verify that the results did not significantly 
change when the geometry was modified from the full fer- 
ric wheel with six iron atoms to the cluster with two iron 
atoms. This was first done with basis set A, so that an 
identical basis set was applied for the full molecule and 
the fragment. 

With the code MOLPRO, calculations at the level of 
MCSCF and MRPT2 were performed. For these meth- 



TABLE I: Geometrical parameters of the molecular cluster, 
in A. The mirror plane is the yz plane. 



Fe 


-1.568236 


2.716264 


0.000000 


0(1) 


0.000000 


1.952334 


-1.010690 


0(2) 


0.000000 


3.375149 


1.066879 


H 


0.000000 


2.019278 


-1.966749 


H 


0.000000 


4.227152 


1.505766 


O 


-1.660933 


4.316006 


-1.146858 


O 


-2.906901 


3.596739 


1.146858 


C(3) 


-2.381312 


5.371681 


-1.023420 


C(l) 


-3.232182 


5.601224 


0.000000 


C(2) 


-3.460321 


4.748871 


1.023420 


H 


-4.021835 


4.973378 


1.723433 


H 


-3.551594 


6.475448 


0.000000 


H 


-2.293279 


5.970805 


-1.723433 


O 


-2.815456 


1.465882 


-0.940540 


O 


-1.442379 


1.009753 


1.120700 


H 


-3.525746 


1.728168 


-1.530737 


H 


-1.465828 


1.108581 


2.107450 


+1 


-2.774296 


-0.177610 


0.182920 



ods, a modified basis set (from now on labeled as ba- 
sis set B) was used. For iron, a [8s5pSd] |3£j basis set 
was chosen, where a /-exponent of 2.48 was added which 
was optimized in a preceding calculation. For oxygen, a 
[4s3p] J33 basis set was chosen where a d-exponent of 0.8 
was added. The basis sets for carbon and hydrogen were 
chosen accordingly to basis set A. Thus, the final basis 
set B was of the size [8s5p3dlf] (iron), [5s^p] (carbon), 
[4s3pld] (oxygen) and [2s] (hydrogen). The enlargement 
of the basis set for iron and oxygen in basis set B com- 
pared to basis set A is a necessary procedure to properly 
account for the needs of post Hartree-Fock calculations, 
which in contrast to the calculations at the UHF and 
B3LYP level include electronic excitations and thus re- 
quire a larger virtual orbital space. To investigate the 
impact of the enlargement of the basis set on the results, 
calculations were carried out with MOLPRO using either 
basis set A or basis set B at the MCSCF level. Despite 
the slightly different basis sets, no significant changes in 
the results were observable, and basis set B can be con- 
sidered as an extension of basis set A for wave function 
based correlation calculations. Subseqently all MRPT2 
calculations were performed with the enlarged basis set 
B. 

The properties of the ferromagnetic (FM) state (all 
spins parallel, total spin 10 /is) and of the antiferromag- 
netic (AF) state (spins alternating up and down, total 
spin 0) were computed with each method. To obtain the 
net charge, a Mulliken population analysis was performed 
at all levels of theory. 



FIG. 2: Spin densities of the molecular cluster for the anti- 
ferromagnetic (AF) state at the UHF level (upper panel) and 
the B3LYP level (lower panel). Both graphs show the spin 
density in the plane given by the planar arrangement of the 
six iron atoms of the primal ferric wheel. For all figures, the 
contour lines range from -0.0004 to 0.0005 in steps of 0.000035 
electrons/(a.u.) 3 . Full lines indicate positive spin density and 
dashed lines indicate negative spin density. 




Depending on the ansatz for the ground state wave 
functions the magnetic coupling in the molecular cluster 
is derived from fits to an Ising or Heisenberg model. 

The analysis in terms of a Heisenberg model 
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Si • S2 — 



(Si+S 2 ) 2 -S 2 -S 2 



is appropriate when the approximate ground state is an 
eigenfunction of (Si + S 2 ) = S 2 . This is the case for MC- 
SCF and MRPT2. The exchange parameter is obtained 
from the difference between the ferromagnetic (S = 5) 
and the antiferromagnetic alignments which is given by 
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FM 
H.tot 
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AF 
H.tot 



-J(Si ■ 
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Si-S 2 ^) 



for quantum spins. The trial wave functions used for 
UHF and B3LYP calculations usually lack invariance un- 
der spin rotation. Since they are constructed as eigen- 
functions of the spin projection S z the energy gain from 
magnetic correlations is analyzed in terms of the Ising 
model: 



Hi 



-J • S\ z ■ S*22 



FIG. 3: Magnetic exchange parameter of the molecular clus- 
ter for different values of the level shift. The upper graph 
(squares) represents the uncorrected values, the lower graph 
(triangles) the corrected values (correction performed accord- 
ing to Roos et al [3(J). 
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III. RESULTS 

Table ITTl summarizes the results for the total energies 
of the ferromagnetic and the antiferromagnetic ground 
states, their differences and the magnetic exchange pa- 
rameter J for the molecular cluster. 



The corresponding energy difference between ferro- 
magnetic (S\ z = S2z) and antiferromagnetic (Sx z — 
—S2z) alignment is given by 



P FM P AF _ 
a I,tot ~ ^I,tot — 



-2-J-Sf z = -12.5- J 



which agrees with the Heisenberg value in the classical 
limit Si — ► 00. 

The experimental value for the exchange parameter of 
the primal ferric wheel was found to be J=-21 K 16] . 
Thus for all results in this dimension the difference in the 
exchange parameter between the Ising and the Heisen- 
berg model is about a few Kelvin; with respect to this 
actuality we consider the Ising model approach as a valid 
description of the magnetic coupling in the cluster. 

We first compare the results of one-determinantal 
methods for the model cluster and the full molecule. In 
a second step, we focus on the influence of correlations. 
In our previous studies of the full molecule [3 , we deter- 
mined an exchange parameter of J= +7 K at the UHF 
level and J=-31 K at the B3LYP level. 

The UHF result for the full molecule is perfectly repro- 
duced by the cluster yielding J=+7 K. A good agreement 
between cluster and periodic system at the Hartree-Fock 



level was already observed in earlier studies, e.g. NiO 15 1, 

23. 



KN1F3 and K2N1F4 Um or Ca 2 Cu0 3 and Sr 2 Cu0 3 
At the B3LYP level, the computed exchange parameter 
is J=-20 K. The agreement with the value for the full 
molecule is thus not as perfect as at the UHF level, but 



still reasonable (a similar deviation was found, for ex- 
ample, when comparing B3LYP exchange couplings for 
K2CUF4 from cluster |4£| and periodic systems [ll|; it 
should be mentioned that there seem to be fewer com- 
parisons between cluster and periodic calculations at the 
B3LYP level) . In addition, the Mulliken population anal- 
ysis fTables ITTT1 and Hvjl demonstrates that the charge is 
practically identical at the various levels: B3LYP popu- 
lations for the cluster and the full molecule [9| are virtu- 
ally identical, and similarly UHF populations agree very 
well. In addition, the MCSCF charge agrees with the 
UHF charge. As a whole, we feel that the cluster model 
can be considered to be physically valid. 

Concerning the individual charge, we note that as a 
confirmation of the former findings 9j , the net charge of 
iron is between 1.55 (B3LYP) and 2.26 (MCSCF) and 
thus far away from the formal charge of +3 of the primal 
ferric wheel [l||- The charge is thus more delocalized 
at the B3LYP level, which can also be seen in the spin- 
density plot in Figure [21 At the B3LYP level the local 
magnetic moment is distributed over the cluster to a cer- 
tain extent, while the level of delocalization is apparently 
smaller at the UHF level. 



To serve as a starting point for the MRPT2 calcu- 
lations, a MCSCF calculation was performed. The ac- 
tive space is made of the iron d-orbitals. The obtained 
value for the exchange parameter is J=-0.5 K. Com- 
pared to UHF, a MCSCF wave function is a better ap- 
proximation to the ground state than a single determi- 
nant, which is the reason for the slight change of the 
exchange parameter from UHF to MCSCF towards the 
experimental value. A direct comparison between the 
calculations with the code CRYSTAL and MOLPRO is 
only approximatively possible, and only in case of high 
spin. First, the MOLPRO energy must be corrected for 
the fact that MOLPRO does not take the interaction 
of point charges into account when computing the to- 
tal energy: when this is done, then the MCSCF energy 
with MOLPRO is -3507.90171+1/(2*2. 774296/. 5291772) 
Eh =-3507. 80634 Eh- The remaining difference to the en- 
ergy computed with CRYSTAL (-3507.80720 E h ) is be- 
cause the MOLPRO MCSCF wave function has identical 
orbitals for up and down spin, whereas the CRYSTAL 
UHF wave function has not; and because the codes use 
different screening parameters for the selection of the in- 
tegrals. 



r 



Within the MCSCF scheme, the wave function is built 
from all Slater determinants representing charge transfer 
configurations between the orbitals in the active space, 
and the orbitals are optimized. The MRPT2 calculation 
is based on the MCSCF orbitals. The reference configu- 
rations are made of all determinants in the active space 
(i. e. the iron ci-orbitals), and subsequently second-order 
perturbation theory is applied. The level of correlation 
can additionally be varied by keeping different sets of 
core orbitals frozen, as described later on in this section. 

While performing the MRPT2 calculations, intruder 
state problems appeared which are due to a near degener- 
acy of the ground state. One possibility to remedy those 
intruder states is to increase the active space, which is 
definitely not possible for the considered molecular clus- 
ter, or, following the proposal by Roos et al y(|, to imple- 
ment a level shift to the MRPT2 calculations. By means 
of this technique, a level shift parameter is added to the 
zcroth order Hamiltonian to avoid those intruder states. 
Thus, the resulting exchange parameters are influenced 
by the level shift (c.f. upper graph (squares) in Fig. |3J 
corresponding to equation 6 in |30| ) . To approximatively 
correct for this effect, a correction to the second order 
energies can be applied afterwards, as was suggested in 
|30|. equation 7. The corresponding data are displayed 
in the lower graph (triangles) in Fig. [3J The quantitative 
dependence of the exchange parameter on the level shift 
is shown in Table Ivl For level shift values smaller than 
0.17 E h {E h = hartree), the MRPT2 results became un- 
stable. The level shift interval (for the corrected values 
of the exchange parameter) from 0.20Eh to 0.30E^ leads 



to an exchange parameter of J=15±2 K, which is a rea- 
sonably stable result. For other systems, similar results 
were obtained, see e.g. de Graaf et al 23] and Hozoi et 
al 29]. 

The results for the MRPT2 calculations are given in 
table lVll Considering the magnetic iron ci-orbitals as the 
only orbitals to be correlated led to an exchange param- 
eter of 

J=-6.4 K. The included configurations are thus the 
charge transfer configurations between the occupied iron 
d-orbitals and the excitations to the virtual orbitals (for 
the high-spin state, only excitations to the virtual or- 
bitals are possible). Therefore, accounting for those 
charge transfer configurations explains the change of 
the exchange parameter from the MCSCF level to the 
MRPT2 level (-0.5K (MCSCF)^ -6.4 K (MRPT2)). 
Adding the 2sp-orbitals of one of the bridging oxygens to 
the orbitals to be correlated gave rise to an exchange pa- 
rameter of J=-4.9 K (0(1)) and J=-7.4 K (0(2)). Tak- 
ing both these 2sp-orbitals into account caused only an 
exchange parameter of J=-5.8 K, i. e. the effect of cor- 
relating both orbital groups is approximatively additive. 

Correlating all oxygen atoms results in a value of -13.7 
K. A further increase of the orbitals to be correlated (up 
to a maximum where only the iron Is-, 2sp-, 3sp-, the 
oxygen Is- and the carbon is-orbitals are kept frozen) 
led to a exchange parameter of J=-14.4 K, which is in 
the range of the experimental value. With respect to the 
geometry of the cluster, the influence of those apical lig- 
and groups to the exchange parameter is the dominant 
one, whereas the oxygen bridge atoms have only little 



TABLE II: Total energies (Eh = hartree), the differences in total energies of the ferromagnetic (FM) and the antiferromagnetic 
(AF) state and exchange parameters J. At the MRPT2 level, all orbitals except the core orbitals were correlated and a level 
shift of O.SEh was applied. The core orbitals include the iron Is-, 2sp-, 3sp-, the oxygen Is- and the carbon Js-orbitals. For 
comparison, the results for the full molecule obtained from ref. Q are given in the first section of this table. 

basis set molecule method FM total energy AF total energy difference of J (K) 
(E h ) (E h ) total energy (mE h ) 



A 


full molecule 


[9] 


UHF 


-13295.73287 


-13295.73125 


-1.62 


+7 


A 


full molecule 


[91 


B3LYP 


-13334.50676 


-13334.51409 


7.33 


-31 


A 


cluster 




UHF 


-3507.80720 


-3507.80693 


-0.27 


+7 


A 


cluster 




B3LYP 


-3514.95130 


-3514.95209 


0.79 


-20 


A 


cluster 




MCSCF 


-3507.90171 


-3507.90175 


+0.04 


-1 


B 


cluster 




MCSCF 


-3508.31083 


-3508.31085 


+0.02 


-0.5 


B 


cluster 




MRPT2 


-3511.04152 


-3511.04209 


0.57 


-14.4 



TABLE III: Charge at various sites, in |e|. Note that the charge is virtually identical for the ferromagnetic and the antifer- 
romagnetic state. For comparison, the results for the full molecule obtained from ref. J9j are given in the first section of this 
table. 



basis set 



molecule 



method O (apical) O (bridge) C(l) C(2)/C(3) 



A 
A 


full molecule [9J 
full molecule [9J 


UHF 
B3LYP 


-1.10 
-0.86 


-1.18 
-0.88 


-0.22 
-0.11 


0.82 
0.62 



A 


cluster 


UHF 


-1.03 


-1.13 


-0.17 


0.77 


A 


cluster 


B3LYP 


-0.80 


-0.89 


+0.01 


0.60 


A 


cluster 


MCSCF 


-1.03 


-1.14 


-0.17 


0.77 


B 


cluster 


MCSCF 


-1.15 


-1.22 


-0.13 


0.85 



TABLE IV: Mulliken charge of Fe, in |e|. Note that the charge is virtually identical for the ferromagnetic and the antiferro- 
magnetic state. The /-population of the iron atoms is negligible. The core orbitals include the iron ls-,2sp-, 3sp-, the oxygen 
Is- and the carbon is-orbitals. For comparison, the results for the full molecule obtained from ref. [ij are given in the first 
section of this table. 



basis set 


molecule 




method 


net charge 


s p d 


A 


full molecule 


[9] 


UHF 


2.16 


6.28 12.28 5.29 


A 


full molecule 


[9] 


B3LYP 


1.56 


6.38 12.37 5.69 


A 


cluster 




UHF 


2.13 


6.29 12.30 5.28 


A 


cluster 




B3LYP 


1.55 


6.39 12.40 5.66 


A 


cluster 




MCSCF 


2.14 


6.29 12.30 5.27 


B 


cluster 




MCSCF 


2.26 


6.37 12.07 5.31 



TABLE V: Total energies (Eh = hartree), the differences in total energies for the ferromagnetic (FM) and the antiferromagnetic 
(AF) state and exchange parameters J for different level shifts at the MRPT2 level for the molecular cluster. All orbitals except 
the core orbitals were correlated. The energies were corrected for the level shift. The core orbitals include the iron Is-, 2sp-, 
Ssp-, the oxygen Is- and the carbon is-orbitals. Basis set B was used. 



level shift FM total energy (Eh) AF total energy (Eh) difference of total energy (mEh) J (K) 



0.20 


-3511.05498 


-3511.05563 


0.21 


-3511.05375 


-3511.05439 


0.22 


-3511.05251 


-3511.05314 


0.23 


-3511.05124 


-3511.05185 


0.24 


-3511.04993 


-3511.05054 


0.25 


-3511.04860 


-3511.04921 


0.26 


-3511.04724 


-3511.04784 


0.27 


-3511.04586 


-3511.04644 


0.28 


-3511.04444 


-3511.04502 


0.30 


-3511.04152 


-3511.04209 



0.65 
0.64 
0.63 
0.62 
0.61 
0.60 
0.59 
0.59 
0.58 
0.57 



-16.4 
-16.1 
-15.9 
-15.6 
-15.4 
-15.2 
-15.0 
-14.8 
-14.6 
-14.4 



TABLE VI: Total energies (Eh = hartree), the differences in total energies for the ferromagnetic (FM) and the antiferromagnetic 
(AF) state and exchange parameters J for different levels of correlation at the MRPT2 level for the molecular cluster. A level 
shift of 0.3Eh was applied. The energies were corrected for the level shift. The core orbitals include the iron Is-, 2sp-, Ssp-, 
the oxygen Is- and the carbon is-orbitals. Basis set B was used. 



level of correlation FM total energy AF total energy difference of J (K) 

(E h ) (E h ) total energy (mE h ) 

(I) iron d-orbitals 

(II) (I) + £sp-orbitals of 0(1) 

(III) (I) + £sp-orbitals of 0(2) 

(IV) (I)+&p-orbitals of 0(1),0(2) 

(V) (I) -fall oxygen orbitals 

(VI) all orbitals except core 



3508.43623 


-3508.43649 


0.26 


-6.4 


3508.63134 


-3508.63154 


0.19 


-4.9 


■3508.64415 


-3508.64444 


0.29 


-7.4 


■3508.84371 


-3508.84394 


0.23 


-5.8 


■3510.58726 


-3510.58781 


0.54 


-13.7 


3511.04152 


-3511.04209 


0.57 


-14.4 



impact on the exchange parameter. As a conclusion, 
for this particular ferric wheel under consideration the 
metal <-> ligand charge transfer configurations dominate 
the metal «-> metal charge transfer configurations over 
the bridging oxygen atoms, which is basically the result 
of the ordering of the magnetic orbitals according to the 
Goodenough-Kanamori rules [41|: essentially, the Fe-O- 
Fe angle is nearly right-angled and thus the coupling is 
small. 

In the experiments, when comparing various ferric 
wheels, an approximatively linear relationship between 
the Fe-O-Fe angle and the value of the exchange coupling 
was observed [8j and confirmed |43. This indicates that 
this angle is crucial for the strength of the coupling (as 
long as the Fe-Fe distance is approximatively constant, 
otherwise this distance may also have an impact). This 
is not in contradiction to the findings here: essentially, 
the strength and the nature of the coupling (ferro- or 
antiferromagnetic) is strongly influenced by this angle, 



but still, to compute the interaction properly, the lig- 
ands must be included in the correlation treatment. This 
was demonstrated, for example, in [24j (figure 2): when 
the correlation treatment is not sufficient, the exchange 
couplings come out too small; but still, the dependence 
on the angle is correct. Even more striking were earlier 
calculations where the ligands were crucial to obtain rea- 
sonable values for the exchange couplings, e.g. for KNiF 3 
and K 2 NiF 4 SI or NiO E3- 



IV. CONCLUSION 

Wave function-based correlation methods were applied 
to a molecular cluster derived from the hexanuclear fer- 
ric wheel [LiFe 6 (OCH 3 )i2-(dbm) 6 ]PF 6 Qjj- The valid- 
ity of the molecular cluster containing two iron atoms 
was tested by means of a one-determinantal approach 
with respect to the formerly calculated results for the 



full molecule [9j for the exchange parameter of the pri- 
mal ferric wheel. In addition, at the UHF and B3LYP 
level, the spin densities and the electronic population, at 
the MCSCF level the electronic population were calcu- 
lated. The population analysis supported the validity of 
the cluster model approach. 

The best result for the exchange coupling parame- 
ter J was obtained at the MRPT2 level (J=15±2 K). 
The influence of intruder state problems on the exchange 
parameter was expli citly investigated, and applying the 
level shift technique [3(| was found to lead to stable re- 
sults. MRPT2 thus gives a more controlled approach 
to the importance of electronic correlations for exchange 
couplings, whereas the density functional results depend 
strongly on the functional chosen. Also, the impact of 
certain atom groups on the exchange parameter was de- 



termined at the MRPT2 level. Correlation of the elec- 
trons of the bridging oxygen atoms was of minor impor- 
tance for the coupling strength. A strong enhancement of 
the computed exchange coupling was however observed 
by additionally correlating the electrons of the apical oxy- 
gen atoms. 
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